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SrRuC>3 is a metallic ferromagnet. Its electrical resistivity is reported for temperatures up 
[ to 1000K; its Hall coefficient for temperatures up to 300K; its specific heat for temperatures up 

ON . to 230K. The energy bands have been calculated by self-consistent spin-density functional theory, 

which finds a ferromagnetic ordered moment of 1.45/ib per Ru atom. The measured linear specific 
heat coefficient 7 is 30mJ/mole, which exceeds the theoretical value by a factor of 3.7. A transport 
mean free path at room temperature of « IOA is found. The resistivity increases nearly linearly 
with temperature to 1000K in spite of such a short mean free path that resistivity saturation would 
be expected. The Hall coefficient is small and positive above the Curie temperature, and exhibits 
CO ' both a low-field and a high-field anomalous behavior below the Curie temperature. 
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I. INTRODUCTION 



SrRu03 has a nearly cubic perovskite structure The combination of good chemical stability, good metallic 

conductivity || , and easy epitaxial growth on various perovskite substrates [|]-[| makes it attractive for possible mul- 
tilayer device applications 0. There is some evidence JsJ that above «800K the structure reverts to cubic. Nominally 
single crystals are therefore likely to be randomly twinned arrangements of the three orientations of orthorhombic or- 
dering. Below T C =160K it is ferromagnetic with a moment at T=0K ||, the largest ordered moment known to 
arise from 4d electrons. Assuming a nominal valence Ru 4+ and configuration 4d(t2 g ) , full Hund's rule spin alignment 
(f t| g , i^2g) predicts a moment 1[Ib and insulating ("half-metallic") behavior of the majority spins. The deficiency 
of moment has received three explanations, (1) collective (rather than atomic) band magnetism fl0|| (2) spin canting 
JlT| , and (3) incomplete alignment of magnetization in different domains due to large magnetocrystalline anisotropy 
and random crystallographic domain alignment p2| . In this paper we offer new transport and thermodynamic data, 
and a band structure calculation which agrees with the band magnetism view of the reduced moment. In particular, 
the specific heat shows a very large Fermi level density of states, which band theory assigns to both majority and 
minority spins, with majority spins slightly more numerous. The majority spin electrons have a small Fermi velocity, 
^ . suggesting that most of the transport currents are carried by minority spin electrons. 



II. ELECTRONIC STRUCTURE 



The electronic structure calculations were done using the self-consistent linear muffin-tin orbital (LMTO) method 
p3[ in the local spin density approximation (LSDA). The calculations for the distorted unit cell, containing 20 atoms, 
were converged using 64 k-points in the irreducible (1/8) Brillouin Zone. The density of states (DoS) and Fermi 
velocities were calculated by tetrahedron integration using a finer k-point mesh. The structural parameters were 
taken from ref. |Q . 

A first calculation was done for the undistorted cubic structure with 5 atoms per cell. The shape of the DoS 
is similar to that for the real unit cell, except for details. The DoS value per atom at the Fermi energy (Ep), is 
somewhat smaller than for the large cell, making the Stoner factor (1.04) and magnetic moment (0.6 /j,b) smaller. 
The total energy is lower by about 25 mRy per formula unit (fu) in the ferromagnetic configuration compared to the 
paramagnetic case. For the remaining discussion we give only the results from the calculations for the large distorted 
cell. 
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Fig. 1 shows the DoS for the two spins. The bands extending from about -8eV to about -3eV relative to Ep are 
primarily of oxygen p parentage. There is a gap of about leV separating these from bands of primarily (74%) Ru d 
parentage. The Fermi level lies about 2eV above the bottom of these bands, filling to about two-thirds occupancy a 
region of high density of states, presumably associated with the t2 g orbitals which are mixed mostly through 7r-overlap 
with oxygen p states. The positions of these levels accord well with the assignments made by Cox et al. Jig ] based on 
UPS measurements. Above the Fermi level comes an increasing amount of Ru s and Sr s and d levels in addition to 
the Ru d states of e g -character. In the paramagnetic case when the exchange splitting is zero, Ep falls near the peak 
in the DoS so that the Stoner factor becomes quite large, 1.39. Thus a ferromagnetic transition is expected, mainly 
due to the Ru-d band, and the magnetic moment is 1.45 [Ib per Ru atom in the spin polarized calculation. The other 
atoms have only small moments, about 0.06 and 0.14 /zb per Sr and O atom respectively. The total energy difference 
between the ferro- and paramagnetic configuration is about 24 mRy per fu. Our calculated moment exceeds the 
zero-field measured value « 1.1/iB P|< There is some uncertainty in the true value of the saturation magnetization. 
Kanbayasi p2| extrapolates a high-field value as great as 2/ie. It is possible that a large magnetocrystalline anisotropy 
P, ^2|Jlq | pins the magnetization perpendicular to the orthorhomic c-axis of the distorted structure, which is randomly 
aligned in the nominally single-crystal sample. 

It is interesting that our calculated moment is so sensitive to the small distortion of crystal structure from cubic. 
This accords with the discovery by Kirillov et al. 17 that optic phonon modes are affected by the magnetic order. 

Table I shows the calculated value of the density of states at the Fermi surface for each atom species and spin 
type, and also the calculated values of the Fermi velocity and Fermi surface area. Three bands of majority spins cross 
the Fermi level, with quite small Fermi velocities, whereas 8 bands of minority spin cross the Fermi level with more 
ordinary size Fermi velocities. The DoS at Ep, N(Ep), is sensitive to changes in the magnetic moment, because of 
the position of Ep on the downslope of majority and upslope on the minority spins as shown in fig. 1. At convergence 
we find the total N(Ep)— 190 states/ (Ry-cell), where 60 % is majority spin. This DoS is concentrated on the Ru 
sites, making the value very high per Ru atom. The DoS is even higher in the paramagnetic state, where the Fermi 
surface lies in a region of very small Fermi velocity. Ferromagnetic ordering diminishes the area of the Fermi surface 
to 75% of the paramagnetic state value, mostly by diminishing the area of the up spin surface. The DoS is reduced to 
58% of the paramagnetic value. The up spin contribution is diminished because of the reduced area of Fermi surface, 
and the down spin contribution is reduced because of the increased Fermi velocity. 

Recently we learned of an unpublished LAPW calculation by Singh [jl8| . Near the Fermi level, Singh's results agree 
very closely with ours, including the size of the ordered moment. Farther away there are interesting differences of 
detail, but reasonable qualitative accord. 



III. SPECIFIC HEAT 



Polycrystalline samples of mass ~ 0.5g were prepared for specific heat measurements. Stoichiometric mixtures of 
SrC>2 and Ru02 (99.9% Atomergic Chemetals) were first heated for 24 h at 1000°C and retired again for 80 h at 
1100°C. Several grindings were done in an agate mortar between the thermal treatments in order to get homogeneous 
samples. After 50 hours at 1150°C the reactions were fully realized, and slow quenching (10°C/h) was done under 
oxygen. X-ray powder diffraction of this product showed a single phase consistent with the previous report of Jones 
et al. f|. 

The low temperature heat capacity measurement used a thermal relaxation technique with 1% accuracy ]l9| . The 
polycrystalline sample, 47mg in mass, was measured in zero external field between 1.2 and 32K. Fig. 2 shows the 
data below 7K in the usual C/T vs. T 2 plot. Below 10K the data are fit by the low T approximation 

^ = 1+ ll R ^Tl (!) 
T 5 © D w 

where 7 is the Sommerfeld constant and ©d the Debye temperature in the limit T — > 0. From a least squares fit 
between 1 and 7K we extract the values 7=6.0mJ/atg (1 atg is 47. 3g or 1/5 mole SrRuOs) and 8d=368K. Deviations 
occur at higher temperature, indicating that Od has a minimum near 25K and then increases with T. Such a 
variation is quite common and reflects additional modes in the lOmeV range beyond those in the Debye distribution. 
The measured 7 corresponds to a Fermi Level DoS of 173 states/(Ry formula unit) (both spins) which exceeds the 
band theoretical value by a factor 1 + A 7 = r y ex p/')th = 3.7. Almost an identical enhancement of experiment over 
theory was recently found for the layered metal Sr2Ru04 |2^,^l|. By contrast, the rutile structure metal Ru02 has a 
theoretical density of states at the Fermi surface which is only half as large per ruthenium atom p2[| . Comparing with 
the measured specific heat [^3j this corresponds to a value A 7 = 0.45 close to what is expected from normal electron- 
phonon mass enhancement. Thus Ru02 appears to be a conventional metal, while both SrRuOs and Sr2Ru04 have 
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large mass enhancements A 7 presumably caused less by electron-phonon effects than by spin fluctuations or some 
other kind of Coulomb correlation effect. 

Fig. 3 shows the high temperature specific heat of a 0.36g sample cut from the same ceramic mass, and measured 
in an adiabatic, continuous heating calorimeter [fj4| . The accuracy is 1% and the scatter 0.02%. The full width of the 
ferromagnetic transition exceeds the temperature steps used in data reduction by about two orders of magnitude. 

The total specific heat approaches 20J/(atg K) at T=250K, a relatively low value compared to the Dulong-Petit 
saturation at 24.94 J/(atg K). Neglecting electronic contributions, this sets a lower limit to the Debye temperature, 
Od(250K)> 540K, significantly higher than the low temperature value, 368K. This comparison suggests that the 
phonon density of states is widely dispersed, including high frequency optical modes from ruthenium-oxygen bond 
stretching vibrations. 

The ferromagnetic transition centered at T c ps 160K appears as a surprisingly small step rather than the quasi- 
logarithmic divergence frequently seen. A closer look at the derivative of the specific heat in fig. 3 shows two transition 
components at 158 and 162K, perhaps caused by sample inhomogeneity. No significant anomaly is seen in the range 
35-45K where the resistivity shows some structure (see below.) 

The entropy associated with the transition can be estimated by subtracting phonon and spinless electronic contri- 
butions to the specific heat. Most of the magnetic entropy in model systems with well- localized spins is released in 
a specific heat peak extending from «0.5T C to «1.5T C , i.e. from 80 to 240K, and should equal i?ln3 or i?ln2 for 
full disordering of localized spin of 1 or 1/2 ions which are fully ordered at T = 0. Since we have no clear way to 
make this subtraction, we used two different methods modeled after a critical transition and a mean-field transition 
respectively. Phonons were modeled by sums of Einstein functions and electrons were assumed to give a contribution 
-fT. The results for the integrated entropy AS* in the interval 40 to 200K were « 1.1 J/(atg K) and 2.7 J/(atg K) 
respectively by the two methods. These correspond to 12% and 30% of i?ln3. Apparently no more than 1/3 of the 
entropy of S — 1 Ru spins is removed by ferromagnetic ordering. However, if we assume an electronic term 7T with 
7 independent of temperature and having the value measured at low T, as was done in both entropy estimates, then 
there is a lot of electronic entropy (« .7R per Ru atom) removed by cooling from 200K to T = Q. It is probably not 
correct to attempt a complete separation between electronic and magnetic entropy. 

IV. TRANSPORT PROPERTIES 

Single crystal samples of dimensions « 1 x 1 x 0.5mm were grown for resistivity measurements by slow cooling 
in SrCl2 flux |25|]. First polycrystalline SrRuC>3 was prepared from stoichiometric quantities of RuC>2 (99.9%) and 
SrC03 (99.99%). The single-phase product was mixed with SrCl2 anhydrous in the weight ratio 1:30, and placed in 
a 50cm 3 platinum crucible covered with a platinum lid in a vertical furnace with programmable temperature control. 
The refractory crucible supports induced a small temperature gradient towards the base of the crucible. The flux 
mixture was heated to 1260°C (200°C/h) for homogenization of the melt, soaked for 40h, and cooled to 800°C (at 
2°C/h). The crucibles were then removed and cooled quickly to room temperature. Crystals a few mm in size were 
easily separated from the frozen flux by immersing the crucible in hot water. 

Resistivity was measured in a four probe configuration on samples of typical dimensions « 1 x 0.5 x 0.05mm 3 . 
Approximately 10 samples were measured. The major features were reproducible, i. e. the break in slope at T c , the 
qualitative shape, and the absolute magnitude (to within a factor of 2). We attempted to measure the anisotropy of 
resistivity between the long axis and the thin axis of the samples; the effect was less than 10%. Given that LSDA band 
theory (section II) finds an orthorhombic anisotropy varying from 10% to 40%, the lack of measurable anisotropy of 
our samples may mean that there is a random distribution of orthorhombic domains or twins in our samples. Results 
for one of our samples up to 1000K are shown in fig. 4. 

The electrical resistivity has been reported previously at temperatures below 300K. Polycrystalline samples were 
reported in refs. P, fl4|]2lj |, single crystals in ref. psj] , and oriented thin film samples in refs. p-fl. Our data agree with 
all measurements concerning the shape of p(T) and the break in slope at T c w 160K. Measurements on polycrystalline 
samples tend to yield absolute magnitudes of p(T) about 5 times larger than seen on single crystal or oriented film 
samples. Our measurements, as well as previous ones on single crystals or oriented films show a surprising sample- 
dependence in the low temperature power law, with the exponent p (defined by p(T) — p(0) cx T p ) appearing to vary 
between 1 and 2. 

Fig. 5 shows the resistivity on a different sample in the range to 300K, along with the derivative dp/ dT and a 
Bloch-Griineisen curve. As is typical in metallic ferromagnets, there is a peak in dp/ dT at the Curie point, similar 
to the peak in C{T). We presume that the resistivity can be composed of a sum of electron-phonon scattering and 
spin-fluctuation scattering. Above the Curie temperature, the spin-fluctuation part should slowly saturate at an 
approximately temperature-independent value. The Bloch-Griineisen curve in fig. 5 is chosen to match the slope of 
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the data at high temperature. Insofar as this curve represents a fair estimate of the electron-phonon contribution 
to p(T), we can use it to extract a value of the electron-phonon coupling constant At r - Lacking more detailed 
information, we assume that the scattering rate is the same for up and down spin electrons, and roughly constant for 
all electrons at the Fermi surface. Then a = ti^r/An where the Drude plasma frequency fi^ is defined by 

^lx = 77— y^vleJiekv - cf). (2) 

"coll , 

KG 

From the values given in table I, the xx, yy, and zz components of the Drude plasma frequency tensor have the 
values 3.0, 2.5, and 2.3eV in the ferromagnetic state, and 2.2, 2.5, and 2.3eV in the paramagnetic state. For a random 
mixture of orthorhombic domains we should use the root mean square value, f2 p =2.6eV (ferromagnetic) or 2.3eV 
(paramagnetic). Using a Debye spectrum with 6_d =368K, the curve of fig. 5 corresponds to A tr =0.5, calculated 
using the paramagnetic state value of the Drude plasma frequency. This is similar to the value found in ref. p2|| for 
R11O2. However, because LSDA band theory agrees well with the specific heat 7 in the latter compound but not in 
SrRuC>3, the value of A tr found here cannot be regarded as particularly reliable. 

It is surprising how rapidly thermal scattering sets in as temperature increases from T=0, much more rapidly than 
the Bloch T 5 -law, but sample-dependent as mentioned above. We do not understand the origin of either the rapid 
rise of p(T) or its sample dependence. Independent of a model for the scattering, however, we can estimate the carrier 
mean free path £, using the equation [^7j 

[A fs ^ + A FSl n] (3) 
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where Ap$ is the Fermi surface area. This result depends only on the assumption of not having too strong a variation 
of £k from point to point on the Fermi surface. For paramagnetic materials with Apst = Apsi, eqn. 3 permits a 
robust estimate of the average mean free path, since LDA band theory generally gets the shape of the Fermi surface to 
very good accuracy, not contaminated by renormalization or other likely errors of LDA. In the ferromagnetic case, we 
are less certain about the correctness of LSD theory concerning spin-splitting; therefore the separate values of Apsi 
and ^4_fsj. are less robust. Using the room temperature resistivity of 200/xJ7cm for the sample shown in fig. 5, and 
the paramagnetic value of the Fermi surface area, the room temperature mean free path is 7 A, as small as is found in 
good superconductors with the A15 structure. However, unlike the case of A15 superconductors, there is no sign of 
"saturation" of the resistivity (p(T) approaching a temperature-independent value.) If the value p(300K) =150/if2cm 
from fig. 4 is used, £ becomes « IOA. This is still short enough that the Boltzmann quasiparticle gas theory should 
not be applicable p7| . 

The origin of this short room temperature mean free path is the small Fermi velocity found by LDA band theory 
in the paramagnetic phase. One way of having a larger mean free path is to argue that the paramagnetic state 
has fluctuating spins with medium range order, so that locally for time intervals not too long the energy bands are 
really the ferromagnetic bands. The areas of the up and down spin Fermi surfaces are not very different. However, 
assuming a roughly constant time between scattering events, the small velocity of the majority spin carriers suggests 
a significantly smaller mean free path for the up spins than for the down spins. Therefore, we assume that the second 
term of eq. 3 dominates, and find from the data of fig. 5 an upper limit ^(300K)=17A, or 22A using the data of fig. 
4. The up spin electrons are unlikely to have a mean free path larger than w 9A By T = 500K, even the down spin 
electrons should have I < 10A, and Boltzmann theory should fail visibly. However, the resistivity as shown in fig. 4 
continues to rise almost linearly for T up to 1000K. The origin of this behavior is not known, but is reminiscent of 
the resistivity seen in high T c superconductors and seen also in the high T metallic phase of VO2 [§8| . 

Our samples show a small but significant enhancement in dp/dT in the range 35-45K (fig. 5). This effect is 
sample-dependent. In exactly this range, Kanbayasi [^9| found a dramatic change in magnetocrystalline anisotropy 
and presumed change in orientation of the magnetization in certain samples called "tetragonal" rather than "pseudo- 
cubic" (orthorhombic). However, no sign of anomaly can be seen in the specific heat measured on polycrystallinc 
samples in this range (fig. 3). Possibly our nominally single crystal samples actually contain varying amounts of 
Kanbayasi's second crystal structure. 

The Hall effect was measured in a six-probe configuration, with a typical current of 10mA and magnetic field up 
to 12T. The sign of the Hall voltage was calibrated with a sample of Bi 2 Sr 2 CaCu208 superconductor. Fig. 6 shows 
typical curves of Hall voltage versus field at temperatures above and below the Curie temperature. Above T c at 190K 
there is a normal weak positive Hall coefficient, w 2x 10~ 10 m 3 /C, which diminishes gradually at higher T. Just below 
T c at 140K, the data could be interpreted as showing the sum of an anomalous Hall effect of positive sign and an 
ordinary Hall effect of negative sign. However, the crossover field (~ 4T) is much too large. The data at 100K show 
clearly that there are three regimes of field, with an anomalous Hall effect occuring in very weak fields, not resolved 
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in our experiment, and then two higher field regimes, each characterized by approximately linear behavior of Vh with 
applied field, but with a slope change near 6T. 

Since the Fermi surface contains both hole- like and electron-like regions, it is very difficult to predict from the 
band structure how Rh should behave. The temperature variation can be rationalized in a model where electron- 
phonon scattering is dominant on the hole-like parts and spin-fluctuation scattering dominates on electron-like parts. 
As temperature increases for T>T C , electron-phonon scattering increases while spin fluctuation scattering saturates. 
This drives conductivity increasingly onto the electron-like parts of Fermi surface. 

ACKNOWLEDGMENTS 

We thank C. H. Ahn, P. Almcras, L. Mieville, J.-M. Triscone, Z.-G. Ye, and K. Yvon for helpful discussions. PBA 
was supported in part by U. S. NSF grant no. DMR-94 17755. 



permanent address: Dept. of Physics, SUNY, Stony Brook, NY 11794-3800, U.S.A. 
[1] A. Callaghan, C. W. Moeller, and R. Ward, Inorganic Chemistry 5, 1572 (1966). 
[2] C. W. Jones, P. D. Battle, P. Lightfoot, and W. T. A. Harrison, Acta Cryst. C45, 365 (1989). 
[3] Y. Noro and S. Miyahara, J. Phys. Soc. Jpn. 27, 518 (1969). 

[4] C. B. Eom, R. J. Cava, R. M. Fleming, J. M. Phillips, R. B. van Dover, J. H. Marshall, J. W. P. Hsu, J. J. Krajewski, and 

W. F. Peck, Jr., Science 258, 1766 (1992). 
[5] C. B. Eom, R. B. Van Dover, J. M. Phillips, D. J. Werder, J. H. Marshall, C. H. Chen, R. J. Cava, R. M. Fleming, and 

D. K. Fork, Appl. Phys. Lett. 63, 2570 (1993). 
[6] X. D. Wu, S. R. Foltyn, R. C. Dye, Y. Coulter, and R. E. Muenchhausen, Appl. Phys. Lett. 62, 2434 (1993). 
[7] L. Antognazza, K. Char, T. H. Geballe, L. L. H. King, and A. W. Sleight, Appl. Phys. Lett. 63, 7 (1993). 
[8] G. L. Catchen, T. M. Rearick, and D. G. Schlom, Phys. Rev. B49, 318 (1994). 
[9] A. Kanbayasi, J. Phys. Soc. Jpn. 41, 1876 (1976). 
[10] J. M. Longo, P. M. Raccah, and J. B. Goodenough, J. Appl Phys. 39, 1327 (1968). 

[11] T. C. Gibb, R. Greatrex, N. N. Greenwood, D. C. Puxley, and K. G. Snowdon, J. Solid State Chem. 11, 17 (1974). 
[12] A. Kanbayasi, J. Phys. Soc. Jpn. 44, 89 (1978). 

[13] O.K Andersen, Phys. Rev. B12 3060 (1975); T. Jarlborg and G. Arbman, J. Phys. F7 1635 (1977). 
[14] M. Shikano, T. K. Huang, Y. Inaguma, M. Itoh, and T. Nakamura, Solid State Commun. 90, 115 (1994). 
[15] P. A. Cox, R. G. Egdell, J. B. Goodenough, A. Hamnett, and C. C. Naish, J. Phys. C: Solid State Phys. 16, 6221 (1983). 
[16] L. Klein, J. S. Dodge, T. H. Geballe, A. Kapitulnik, A. F. Marshall, L. Antognazza, and K. Char, Appl. Phys. Lett. 66, 
2427 (1995). 

[17] D. Kirillov, Y. Suzuki, L. Antognazza, K. Char, I. Bozovic, and T. H. Geballe, Phys. Rev. B51, 12825 (1995). 
[18] D. J. Singh, (unpublished). 

[19] D. Sanchez, A. Junod, J.-Y. Genoud, T. Graf, and J. Muller, Physica C200, 1 (1992). 

[20] Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. Fujita, J. G. Bednorz, and F. Lichtenberg, Nature 372, 532 (1994). 
[21] T. Oguchi, Phys. Rev. B51, 1385 (1995); D. J. Singh, Phys. Rev. B (in press). 

[22] K. M. Glassford and J. R. Chelikowsky, Phys. Rev. B49, 7107 (1994). J. H. Xu, T. Jarlborg, and A. J. Freeman, Phys. 

Rev. B40, 7939 (1989). 
[23] B. C. Passenheim and D. C. McCollum, J. Chem. Phys. 51, 320 (1969). 

[24] A. Junod, E. Bonjour, R. Calemzuk, J.-Y. Henry, J. Muller, G. Triscone, and J. C. Vallier, Physica C211, 304 (1993). 

[25] R. J. Bouchard and J. L. Gillson, Mat. Res. Bull. 7, 873 (1972). 

[26] J. J. Neumeier, A. L. Cornelius, and J. S. Schilling, Physica B198, 324 (1994). 

[27] P. B. Allen, in Quantum Theory of Real Materials, edited by J. Chelikowsky and S. G. Louie (Kluwer, in press.) 
[28] P. B. Allen, R. M. Wentcovitch, W. W. Schulz, and P. C. Canfield, Phys. Rev. B48, 4359 (1993). 
[29] A. Kanbayasi, J. Phys. Soc. Jpn. 41, 1879 (1976). 

FIG. 1. Electronic density-of-states of ferromagnetic SrRu03. Majority spin is plotted upward, minority spin downward. 
The cell contains four formula units. 



FIG. 2. Low temperature specific heat C/T plotted versus T 2 for SrRuOs 
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FIG. 3. Specific heat C(T) (diamonds) and derivative dC(T)/dT (full line; note the scale factor) for SrRu0 3 in the 15-250K 
range 

FIG. 4. Resistivity versus temperature for SrRuOs 

FIG. 5. Resistivity versus temperature for SrRuOs from to 300K, with a Bloch-Griineisen curve (solid line) for comparison, 
and dp/dT in arbitrary units 

FIG. 6. Hall voltage versus field for SrRuOs at three temperatures 
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TABLE I. Calculated properties of SrRuQ 3 
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